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Abstract 

The Standard Model includes neutrinos as massless particles, but neutrino oscillations 
showed that neutrinos are not massless. A simple extension of adding gauge singlet fermions 
to the particle spectrum allows normal Yukawa mass terms for neutrinos. The smallness of 
the neutrino masses can be well understood within the seesaw mechanism. We analyse two 
cases of the minimal extension of the standard model when one or two right-handed fields 
are added to the three left-handed fields. A second Higgs doublet is included in our model. 
We calculate the one-loop radiative corrections to the mass parameters which produce mass 
terms for the neutral leptons. In both cases we numerically analyse light neutrino masses 
as functions of the heavy neutrinos masses. Parameters of the model are varied to find light 
neutrino masses that are compatible with experimental data of solar Atoq and atmospheric 
Am^^jjj neutrino oscillations for normal and inverted hierarchy. 
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1 The model 



We extend the Standard Model (SM) by adding a second Higgs doublet and right-handed 
neutrino fields. The Yukawa Lagrangian of the leptons is expressed by 
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k=l 



Cy = -Y, WRTk + ^bijAfc) Dl + H.c. (1.1) 



in a vector and matrix notation, where = ^'^2^%- In expression (1.1) the Ir, uji, and 
Dl = {t^l ^l)^ are the vectors of the right-handed charged leptons, of the right-handed neutrino 
singlets, and of the left-handed lepton doublets, respectively, and $q,, a = 1, 2 are the two Higgs 
doublets. The Yukawa coupling matrices are hl x n^, while the are X- ni- 

In this model, spontaneous symmetry breaking of the SM gauge group is achieved by the 
vacuum expectation values ($fe)vac = vjjV^. ^ = 1>2. By unitary rotation of the Higgs 
doubletswe can achieve (^5)vac — ^ 

/^/2 > and ($§)vac = with v ~ 246 GeV. The charged- 
lepton mass matrix Mi and the Dirac neutrino mass matrix M/j are 

Me = ^ri and = Ai , (1.2) 

with assumption that Mi = diag {me,m^, m,-). The mass terms for the neutrinos can be written 
in a compact form with an (n/, -|- Ufj) x (n/, -|- npj) symmetric mass matrix 

where the hat indicates that Mr is a diagonal matrix. Mj^ can be diagonalized as 

U'^M^ U = m = diag (mi, m2, . . . ,m„^+„^) , (1.4) 

where the rrii are real and non-negative. In order to implement the seesaw mechanism [U [2] 
we assume that the elements of are of order mo and those of M/j are of order m/j, 
with rriD mR. Then the neutrino masses rrii with i = 1,2,. ..,ni are of order m'jj/mji, 
while those with i = til + 1, . . . ,nL + hr are of order rriR. It is useful to decompose the 
(riL+nR) X {riL+nR) unitary matrix U asU = ^| ^ , where the submatrix Ul is x (riL+nR) 
and the submatrix Ur is ur x (hl + ur) [3l Uj. With these submatrices, the left- and right- 
handed neutrinos are written as linear superpositions of the ni+nR physical Majorana neutrino 
fields Xi'- = UlPlx and vr = UrPrx, where Pl and Pr are the projectors of chirality. 

We can express the couplings of the model in terms of mass eigenfields, where three neutral 
particles are coupling to neutrinos. Interaction of the Z boson with the neutrinos is given by 



£M = ^ Z,xr [Pl (uIUl) - Pr (C/Jt/2)J x , (1-5) 
where g is the SU(2) gauge coupling constant and c«, is the cosine of the Weinberg angle. 
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The full formalism for the scalar sector of the multi-Higgs-doublet SM is given in Refs. ^i4j. 
The Yukawa couplings of the Higgs bosons Hjj to the neutrinos are given by 



(1.6) 

with A;, = 6fcAfe, where b are two-dimensional complex unit vectors. The neutral Goldstone 
boson G^^ is given by the vector bz with bz = («,0). 

Once the one-loop corrections are taken into account the neutral fermion mass matrix is 
given by 

mW=( ME + 6Ml\/ 6Ml ^^^^^ 



Md + SMd Mr + 6 Mr J \ Md Mr 



where the Osxs matrix appearing at tree level (1.3) is replaced by the contribution 6Ml. This 



correction is a symmetric matrix, it dominates among all the sub-matrices of corrections. Neg- 



lecting the sub-dominant pieces in (1.7), one- loop corrections to the neutrino masses originate 



via the self-energy function Sf (0) = ^^(0) + sJ^'^"^(0) + i;f ^°^(0), where the ^'^°'^°^(0) 
functions arise from the self-energy Feynman diagrams and are evaluated at zero external mo- 
mentum squared. In the calculation of the self energies the neutrino couplings are determined 



by eqs. (1.5) and (1.6). Each diagram contains a divergent piece but when summing up the 
three contributions the result turns out to be finite. 

The final expression for one-loop corrections is given by [5] 



where the sum runs over all neutral physical Higgses Hjj. 

2 Case nn = 1 

First we consider the most minimal extension of the standard model by adding only one right- 
handed field I'R to the three left-handed fields contained in ly^. 

We parametrize Ai = ^^^^ and A2 = ^^^^^02^ with |al| = 1 and |al| = 1. The 



symmetric mass matrix Mi, (1.3) in block form is diagonalized by 



U^MM = ( "^f' \V={^'' ° ] • (2.9) 
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The non zero masses in A4i and Adh are determined analytically by finding eigenvalues of the 
hermitian matrix M^mI. These eigenvalues are the squares of the masses of the neutrinos 
Ml = diag(0,0,mi) and = m/^. The solutions m|) = mhmi and m|j = {mh — ~ 
correspond to the seesaw mechanism. 

We can construct the diagonalization matrix U for the tree level from two diagonal matrices 
of phases and three rotation matrices Utrec = f^0(</'j)f^i2(ai)f^23(a2)f^34(/?)f^ii where the angle 
j3 is determined by the masses mi and nrih. The values of and Oj can be chosen to cover 
variations in M^). 

For calculation of radiative corrections we use the following set of orthogonal complex 
vectors: bz = (j,0), hi = (1,0), 62 = (0, i) and 63 = (0,1). Diagonalization of the mass 
matrix after calculation of one-loop corrections is performed with a unitary matrix f/ioop = 
UegvUip{ipi, ip2, ifs), where C/ggv is an eigenmatrix of Alj)^^ m!)^^'^ and is a phase matrix. The 
second light neutrino obtains its mass from radiative corrections. The third light neutrino 
remains massless. 

It is possible to estimate masses of the light neutrinos from experimental data of solar 
and atmospheric neutrino oscillations ^ assuming that the lightest mi^ = 0. Considering the 
normal ordering of the light neutrinos we receive nii-^ = 5.0 it 0.2 x 10^^^ GeV and mi^ = 
8.7 lb 0.3 X 10^^^ GeV. Numerical analysis shows that we can reach those values for a heavy 
singlet with a mass bigger than 830 GeV, see Fig. [T} 
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Figure 1 : Calculated masses of two light neutrinos as a function of the heavy neutrino mass ruh . 
The mass of the third light neutrino is zero, when nji = 1. Solid lines show the boundaries of the 
allowed neutrino mass ranges when the model parameters are constrained by the experimental 
data on neutrino oscillations. The purple arrows indicate the values of mi^ which do not satisfy 
the allowed experimental neutrino mass ranges. Due to the scale, the allowed and mi^ 
values form a band, their scattered values are shown separately in the right plots. 
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3 Case nR = 2 



If we add two singlet fields to the three left-handed fields vl^ the radiative corrections give 
masses to all three light neutrinos. 

■ A \/2 I ^Do^l"^ \ A \/2 I ^Do'^^ \ ■ 1 I -> I 

Mow we parametrize Ai = — and A2 = — \ with oi = 1, 

\hi\ = 1, \a2\ = 1 and I62I = 1- Diagonalizing the symmetric mass matrix (1.3) in block 
form we write: 

U-M.U = U-\ m,2^ \u=('i' l )- (3-10) 



The non zero masses in M; and M/i are determined by the seesaw mechanism: m-|,. ~ rrih^mi^ 
and m|j. ~ iT^h ^ « = 1, 2. Here we use the ordering of masses mi > 1712 > m^. The third light 
neutrino is massless at tree level. 

The diagonalization matrix for tree level C/tree = f^i2(oi; «2)f^egv(/3i)C^0(0j) is composed of 
a rotation matrix, an eigenmatrix of C/]^Mi,Mii{7j*2 ^i^d ^ diagonal phase matrix, respectively. 

For calculation of radiative corrections we use the same set of the orthogonal complex vectors 
bi as in the first case. Diagonalization of the mass matrix including the one- loop correction 
is performed with a unitary matrix C/ioop = UcgvU^i^i), where C/egv is the eigenmatrix of 
mI)^^ mI)^^'^ and U,p is a phase matrix. 

In numerical calculations the model parameters as well as the derived masses of the light 
neutrinos are obtained in several steps. First, the diagonal mass matrix for the tree level is 
constructed. The lightest neutrino is massless, and the masses of other two light neutrinos are 
estimated from experimental data on solar and atmospheric neutrino oscillations. The masses 
of the heavy neutrinos are input parameters. This diagonal tree level matrix is used to constrain 
the parameters Oi and (pi that enter the tree-level mass matrix and its diagonalization matrix 
Utree- Then the diagonalization matrix is used to evaluate one- loop corrections to the mass 
matrix. Diagonalization of the corrected mass matrix yields masses for three light neutrinos. If 
the calculated mass difference is compatible with the experimental neutrino mass difference, the 
parameter set is kept. Otherwise, another set of parameters is generated. Figure [2] illustrates 
the obtained results. Both normal and inverted neutrino mass orderings are considered. 
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Figure 2: The masses mi- of the hght neutrinos as functions of the heaviest right-handed 
neutrino mass m/i^, for the case hr = 2. The value of the second heaviest right-handed 
neutrino mass is fixed at = 100 GeV. The plot on the left represents the normal hierarchy, 
and the plot on the right represents the inverted hierarchy of the light neutrinos. The wide solid 
lines indicate the area of the most frequent values of the scatter data. The nearly degenerate 
masses mi^ and rrii^ are shown separately in the lower right plot for the case of the inverted 
hierarchy. 

4 Conclusions 

For the case ur = 1 we can match the differences of the calculated light neutrino masses to 
AuIq and Am^^j^ with the mass of a heavy singlet larger than 830 GeV. Only normal ordering 
of neutrino masses is possible. 

In the case riR = 2 we obtain three non vanishing masses of light neutrinos for normal 
and inverted hierarchies. The numerical analysis shows that the values of the light neutrino 
masses (especially of the lightest mass) depend on the choice of the heavy neutrinos masses. 
The radiative corrections generate the lightest neutrino mass and have a big impact on the 
second lightest neutrino mass. 

In future we plan to apply our parametrization to study the r polarization coming from the 
decay of a ly boson in the data of the CMS experiment at LHC and thus determine restrictions 
to the parameters of the neutrino sector. 
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